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摘 要 : 火星 是 重要 的 地 外 天 体 探 测 目标 之 一 ， 对 火星 表面 进行 的 探测 、 研 究 表 明火 星 表 
面色 经 存在 液态 水 ,水 是 生命 存在 的 基础 与 前 提 条 件 , 因此 在 次 表层 寻找 不 同形 式 的 水 是 目 
前 火星 探测 的 重要 科学 目标 之 一 。 近 十 七 年 来 ,欧洲 火星 快车 上 搭载 的 火星 次 表层 和 电离 层 
探测 先进 雷达 以 及 美国 火星 勘测 轨道 器 上 搭载 的 浅 表 层 雷 达 开 展 了 对 火星 次 表层 的 大 量 探 
JU ,由 于 雷达 的 工作 原理 ， 雷达 获取 的 原始 回 波 数据 需要 经 过 距离 向 、 方 位 向 处 理 、 电 离 层 
校正 等 数据 处 理 步骤 才能 完成 成 像 ， 从 成 像 数据 中 提取 科学 信息 ， 取 得 了 众多 科学 成 果 。 火 
星 次 表层 探测 雷达 的 数据 处 理 技术 在 火星 次 表层 探测 上 发 挥 了 不 可 或 缺 的 作用 ,具有 很 强 的 
代表 性 与 参考 价值 。 综述 目前 利用 雷达 对 火星 次 表层 展开 的 探测 与 研究 ,介绍 了 火星 次 表层 
探测 雷达 数据 处 理 技术 ,列举 了 取得 的 部 分 科学 成 果 ,并 展望 将 要 投入 使 用 的 火星 次 表层 探 
测 雷 达 。 
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0 引言 


火星 是 太阳 系 中 的 类 地 行星 之 一 , 相 比 起 同 为 类 地 行星 的 水 星 、 金 星 , 火星 在 大 气 、 
地 貌 等 方面 与 地 球 相 似 度 较 高 , 是 太阳 系 内 与 地 球 相似 度 最 高 的 行星 。 从 人 类 进入 航天 时 代 
开始 , 火星 就 是 最 重要 的 地 外 天 体 探测 目标 之 一 , 美 苏 以 及 其 他 国家 总 计 进行 了 45 次 火星 飞 
越 、 绕 飞 、 着 陆 以 及 巡视 的 探测 尝试 趾 ， 仅 次 于 "对 月 球 的 124 次 。 火 星 探 测 的 第 一 次 高 峰 期 
出 现在 1960-1975 年 , 在 此 期 间 水 手 4 号 、6 号 、7 号 对 火星 进行 了 飞越 探测 , 近 距 离 拍摄 火星 ， 
为 人 类 提供 了 火星 表面 的 清晰 图 像 。 水 手 9 号 首次 完成 了 对 火星 的 绕 飞 探测 ， 而 海盗 1 号 、2 
号 在 火星 表面 完成 了 软 着 陆 中 ， 同 期 苏联 也 发 射 了 多 个 探测 器 进行 火星 探测 趾 。1992 年 之 后 
火星 探测 再 次 迎 来 高 峰 ， 火 星 快 车 (Mars Express) 、 火 星 勘察 轨道 器 (Mars Reconnaissance 
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27 Orbiter, MRO)、 火 星 全 球 勘 测 者 (Mars Global Survyor MGS) 、 凤 凰 号 着 陆 器 、 勇 气 号 、 
28 机遇 号 、 好 奇 号 火星 车 等 探测 器 都 抵达 火星 开始 探测 任务 ， 对 火星 进行 了 多 波段 光学 遥感 、 
29 ”雷达 次 表层 探测 、 着 陆 巡 视 等 多 种 方式 的 探测 ,绘制 了 火星 表面 的 地 貌 细节 、 矿 物 分 布 、 表 
30 ” 面 亮 温 分 布 ， 对 火星 次 表层 结构 进行 了 研究 ， 在 火星 上 和 寻找 液态 水 及 其 存在 过 的 证 据 。 

31 对 火星 表面 进行 的 一 系列 探测 结果 表明 , 火星 表面 曾经 有 液态 水 存在 , 而 在 火星 南极 冰 
32 ”六 进行 的 次 表层 探测 甚至 为 液态 水 的 存在 提供 了 直接 观测 证 据 ， 欧 空 忆 (ESA) KERER 
33 ” 测 器 上 搭载 的 MARSIS 雷 达 于 2012 年 至 2015 年 间 在 火星 南极 冰 盖 区域 采 集 到 的 数据 显示 , UK 
34 ”六 表 面 下 1.5 千 米 处 存在 宽度 达到 20 千 米 的 液态 水 湖 趾 ,. 火 星 上 存在 不 同 状态 的 水 这 一 发 现 激 
35 ”发 了 人 类 寻找 火星 生命 存在 痕迹 的 热情 。 

36 在 火星 搜索 各 种 形态 水 的 存在 以 及 生命 痕迹 的 过 程 中 , 一 种 先进 的 探测 仪器 发 挥 了 极 大 
37 ”作用 ， 这 种 设备 被 称 为 轨道 器 次 表层 探测 雷达 ， 与 车 载 次 表层 探测 雷达 相 比 ， 星 载 雷 达 最 大 
38 ”的 特点 是 使 用 了 合成 孔径 技术 , 通过 对 图 像 进 行 顺 轨 方向 的 处 理 , 在 离 地 面 数 百 公里 的 高 度 
39 ”上 仍 能 取得 较 高 的 方位 向 分 辩 率 。 轨 道 器 次 表层 探测 雷达 的 探测 深度 可 达 数 百 米 甚至 几 千 米 ， 
日 其 分 辨 率 较 低 , 一 般 为 数 十 米 ,目前 只 有 欧 空 忆 (ESA )2003 年 发 射 的 火星 快车 探测 器 GMars 
41 ”Express》- 和 美国 航空 航天 局 (NASA)〉 2005 年 发 射 的 火星 勘测 轨道 器 上 搭载 了 此 类 仪器 ， 分 
42 别 简称 为 MARSIS 和 SHARAD > MARSIS & f& X Mars Advanced Radar for Subsurface 
43 and Ionosphere Sounding〈 火 星 地 表 和 电离 层 探 测 先 进 雷 达 ) ，SHARAD 全 称 为 Shallow 
44 Subsurface Radar〈 浅 表层 雷达 ) 。 而 在 我 国 “ 剥 娥 ” 探 月 工程 中 大 放 异 彩 的 车 载 次 表层 探 
45 ” 测 雷 达 目 前 还 没有 随 火 星 车 或 者 着 陆 器 登陆 过 火星 ， 开 展 相 应 的 次 表层 探测 。 

46 我 国 的 首次 火星 探测 任务 已 于 2016 年 正式 立项 ，2020 年 7 月 23 日 于 海南 文昌 发 射 场 由 长 
47 ” 征 5 号 运载 火箭 发 射 ， 探测 器 包括 轨道 器 和 巡视 器 四， 轨道 器 与 巡视 器 上 都 搭载 了 次 表层 探 
48 ” 测 雷 达 。 美 国 Mars 2020 任 务 探测 器 于 2020 年 7 月 30 日 由 Atlas V 541 火 第 在 卡 纳 维 拉 尔 角 空 军 
49 ”基地 发 射 升 空 ， 阿 联 疾 的 火星 探测 任务 “希望 号 ”探测 器 于 2020 年 7 月 20 日 由 H-IIA 火 第 在 种 
50 。 子 岛 航天 中 心 发 射 升 宝 ， 欧 洲 与 俄罗斯 的 联合 火星 探测 任务 ExoMars 原 定 于 2020 年 发 射 ， 由 
51 于 疫情 等 多 种 原因 推迟 至 2022 年 ,火星 探测 又 将 迎 来 一 次 高 峰 ， 且 Mars 2020 任 务 和 ExoMars 
52 ”任务 中 近视 器 都 搭载 了 车 载 次 表层 探测 雷达 。 如 果 各 国 的 探测 器 都 能 顺利 抵达 火星 展开 探测 
53 则 如 何 处 理 、 分析 、 解读 火 星 次 表层 探测 雷达 的 数据 将 成 为 火星 次 表层 探测 的 一 大 重点 、 难 
54 Ro 本文 综述 了 国内 外 利用 雷达 探测 技术 进行 火星 次 表层 探测 的 发 展 成 果 与 技术 进展 , 对 于 
55 ”未 来 更 好 利用 火星 探测 任务 提供 的 雷达 数据 有 一 定 意义 。 
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56 
57 1 次 表层 探测 雷达 对 火星 探测 技术 发 展 

58 

59 雷达 发 射 的 电磁 波 穿 透 性 强 ， 且 不 依赖 于 光照 ， 上 共有 全 天 候 性 ,在 月 球 的 次 表层 探测 中 


60 ”发 挥 了 很 大 作用 ，1972 年 搭载 于 阿波 罗 17 号 的 阿波 罗 月 球 探测 仪 、2007 年 日 本 SELENE 拱 载 
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的 月 球 雷 达 探 测 仪 以 及 录 娥 3 号 、4 号 月 球 车 上 搭载 的 测 月 雷达 都 对 月 球 进行 了 次 表层 探测 ， 
并 在 澄海 、 危 海 与 雨 海 等 地 都 探测 到 了 次 表层 分 层 结构 中， 嫦娥 3 号 月 球 车 “ 玉 免 ”与 妙 娥 4 
号 月 球 车 “ 玉 免 2 号 ”都 搭载 了 配置 相同 的 月 球 穿 透 雷达 (Lunar Penetrating Radar, LPR) , 


工作 频段 为 40-80 MHz 与 230-750 MHz， 最 大 探测 深度 超过 100 米 ， 并 拥有 优 于 1 米 的 距离 分 
辩 率 ， 进 行 巡 视 路 线 上 月 球 次 表层 与 月 壤 结 构 探 测 趾 。3 号 雷达 在 月 球 雨 海 的 巡视 区 中 展开 


探测 ， 根 据 探测 结果 估计 风化 层 深度 为 4-6 米 ， 并 确定 约 330 米 深 处 存在 月 岩层 外 。 嫦娥 4 号 
搭载 的 雷达 在 月 球 背 面 的 冯 卡 门 撞击 坑 获 得 了 远 深 于 3 号 着 陆 区 的 探测 深度 ， 并 探 明 了 近视 


区 的 三 个 次 表层 分 层 结构 : 风化 层 〈 深 度 12 米 ) 、 骨 入 岩石 的 粗糙 材料 《深度 12-24 米 ) 以 
及 粗糙 与 精细 物质 的 交替 层 〈 深 度 24-40 米 ) 外。 火星 探测 距离 远 ， 能 量 需求 大 ， 难 度 较 高 ， 


早期 的 火星 探测 一 直 以 光学 为 主 , 雷达 较 少 被 用 于 火星 探测 , 目前 地 


雷达 观测 和 轨道 器 次 


基 
表层 探测 雷达 探测 已 被 正式 用 于 火星 探测 , 已 登陆 火星 的 巡视 器 上 都 没有 搭载 次 表层 穿 透 雷 


达 ， 我 国 的 火星 车 次 表层 探测 雷达 与 美国 Mars 202045 MES “BLAS” (Perseverance) 
火星 车 搭载 的 RIMFAX 中 将 诞生 第 一 台 登 陆 火 星 的 车 载 次 表层 探测 雷达 。 

地 基 雷 达 探 测 是 雷达 技术 应 用 于 火星 探测 的 最 初 方式 , 进行 观测 的 设备 为 位 于 美国 波多 
黎 各 的 阿 雷 西 博 CArecibo) 305 米 望远镜 以 及 位 于 美国 加 利 福利 亚 的 金石 《goldstone) 70 米 


望远镜 nan。 地 


射 


基 雷 达观 测 基于 随机 长 编码 延迟 -多 普 勒 技术 以 及 干涉 延迟 -多 普 勒 技术 


OO, 在 火星 冲 日 时 对 火星 表面 进行 观测 ， 获 取 了 火星 表面 的 高 程 、 反 射 率 、 粗糙 度 等 信息 ， 


并 进行 表面 物质 的 介 电 常 数 反 演 。 


由 于 火星 的 快速 自转 , 入 射 的 电磁 波 会 因 多 普 勒 频 移 发 生 严 重 的 频率 偏 移 ， 大 大 提高 探 
测 难度 "站 ， 且 地 基 雷 达 探 测 技术 无 法 进行 次 表层 探测 。 在 火星 探测 进入 绕 飞 与 着 陆 探测 阶 


段 后 ， 地 基 和 雷达 探测 


更 多 发 挥 着 陆 区 选取 等 辅助 功能 。 


动 了 次 表层 探测 雷达 


在 火星 探测 确立 了 在 火星 次 表层 寻找 不 同形 态 的 水 的 科学 目标 后 ,次 表层 探测 的 需要 推 


技术 的 发 展 , 并 逐渐 在 火星 探测 器 的 载荷 中 占 到 越 来 越 重要 的 地 位 。 目 


前 ,火星 次 表层 探测 雷达 已 获取 了 基本 歼 盖 火星 全 球 的 探测 数据 , 并 识别 出 一 些 水 冰 存 在 的 


区 域 。 


1.1 轨道 器 次 表层 探测 雷达 


次 表层 探测 雷达 拥有 穿 透 火星 地 表 ， 深 入 火星 次 表层 探测 其 结构 的 能 力 。 雷 达 发 射 的 


电磁 波 能 够 在 一 定 程 


BE E36 EHE. UK. AAEN 并 在 两 种 介 电 常数 不 同 的 介质 界面 上 


发 生 反 射 、 散 射 和 透射 ， 这 种 特点 可 以 帮助 雷达 控 明 探测 区 域 中 的 地 质 分 层 ， 比 如 干燥 的 土 


壤 与 冰 层 的 界面 , 冰川 与 基 岩 的 分 层 , 孔 隐 中 充满 冰 或 液态 水 的 多 孔 岩 石 以 及 冰 盖 下 液态 水 


的 存在 等 。 


IL 


已 磁 波 在 各 种 介质 中 并 非 无 损耗 传输 ,电磁 波 在 穿越 介质 时 会 由 于 介质 吸收 而 发 


生 衰减 , 一 种 介质 对 | 


人 磁 波 的 吸收 能 力 一 般 与 这 种 介质 的 介 电 常数 和 电 做 波 的 波长 有 关 , 波 


长 越 长 ， 频 率 越 低 的 


电磁 波 在 同 种 介质 的 中 传播 ， 衰减 的 程度 越 小 ， 穿 透 能 力也 越 强 ， 能 探 


比 ， 次 表层 探测 雷达 工作 频率 一 般 在 1 MHz-1 GHz 之 间 选 取 ， 而 


雷达 的 带宽 决定 了 分 辩 率 ， 带 宽 越 小 ,深度 上 的 分 辩 率 越 高 。 为 了 实现 对 次 表层 进行 高 精度 


分 别 为 1.3-2.3 MHz. 2.5-3.5 MHz. 3.5-4.5 MHz、4.5-5.5 MHz, 真 


SHARAD 使 用 的 工作 频段 为 15-25 MHz， 带 宽 10 MHz， 波 长 范 
SHARAD 由 于 工作 参数 的 差异 ， 分 辨 率 、 探 测 深度 等 和 
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95 ” 测 到 更 深 处 的 地 质 结 构 。 因 |] 
96 

97 ”的 雷达 探测 ， 并 满 / 
98 ”表层 探测 模式 下 使 月 
99 
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Parameter 


Orbit altitude (km) 
Antenna length (m) 
Frequency band (MHz) 
Bandwidth (MHz) 
Transmitter Power (W) 
Pulse length Cus) 
Pulse repetition frequency (Hz) 
Range resolution (m) 
Horizontal resolution in along-track(km) (km) 
Horizontal resolution in cross-track(km) (km) 


Penetration depth (km) 
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MARSIS 


250-900 


40 


1.3-5.5 


127 


150 (vacuum) 


5-10 


107 1.2 轨道 器 次 表层 探测 雷达 原理 
108 在 次 表层 探测 中 ， 电 磁 波 从 环绕 器 


109 ” 层 会 对 电磁 波 的 幅度 与 相位 产生 
反射 、 散 射 ， 最 终 被 雷达 天 线 接 收 ， 还 有 一 部 分 穿 透 火 旺 
外 常数 不 同 的 次 表层 分 层 界面 处 发 生 反射 。 假 设 轨 道 器 轨 证 
112 ” 则 在 发 射电 磁 波 后 时 延 to=2H/e 时 将 接收 到 来 自 星 下 点 火星 
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MARSIS 与 SHARAD 关键 性 能 参数 00.19 


上 安装 的 天 线 发 射 ， 穿 越 火 
站 ， 随 后 电厂 波 抵达 火星 表面 ,一 部 分 电磁 波 被 表面 


E 一 些 区 别 ( 见 表 1) 


探测 深度 的 要 求 ， 雷 达 的 工作 参数 设计 受到 了 很 多 限制 。MARSIS 在 次 
[ 作 频 段 为 1.3-5.5 MHz, 被 划分 为 4 个 子 频 段 , 每 个 频段 带宽 为 1 MHz, 
E 中 波长 范围 从 55-230 m; 
围 为 12-20 m. MARSIS 与 


o 


115 ”to 十 2h/v。 由 于 雷达 的 波束 宽度 较 宽 ， 
116 ” 射 回 波 还 包括 星 下 点 外 表面 散射 杂 波 , 表面 杂 波 存在 一 定时 间 延 迟 , 可 能 与 星 下 点 次 表层 回 


Table.1 Main performance parameters of MARSIS and SHARAD ^"! 


SHARAD 


250-320 


700.28 


15(vacuum) 


0.3-1 


3-6 


1 


的 大 气 层 ， 其 中 电离 


表面 进入 次 表层 ， 穿 透 介 质 ， 在 介 
高 度 为 H，c 为 真空 中 的 光速 ， 
I 回 波 ， 穿 透 地 表 进 入 次 表 
以 小 于 真空 中 的 速度 v=cmm 传 播 (pn 为 介质 的 折射 率 ) ， 假 设 


则 来 自 这 个 界面 的 回 波 抵达 天 线 的 时 间 为 


乾 兰 的 地 表面 积 较 大 ， 表 国 


回 波 除了 星 下 点 表面 反 
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波 信号 重合 , 因此 星 下 点 次 表层 回 波 信 号 经 过 介质 衰减 可 能 被 表面 杂 波 淹没 , 为 了 在 表面 杂 
波 的 存在 下 有 效 识别 次 表层 回流， 必须 对 回 波 进 行 脉冲 压缩 ， 且 学 辩 必须 远 低 于 主因 ,才能 
实现 次 表层 探测 。 


iul 


Surface echo 


一 -ad 一 
| ^» E f \ "d 
6 = A 
E | Sub-surface echo 
H | 
2 | 
N A = C 山 
\ Á 
ANIN VC% v 
Range of compressed echo 
CA) (B) 
© 


(1) Nadir surface echo 

(2) Off-Nadir surface echo 
3 (S) Sub-surface echo 

(4) Off-Nadir surface echo 


T 009 © 
AAA A 
Range 


(C) 
图 1 (A) ABR AGRI TIA RN IO, B) 来 自 表面 与 次 表层 的 回 波 在 时 间 / 距 离 上 的 分 布 示意 
图 ， 强 大 的 表面 回 波 可 能 会 掩盖 较 弱 的 次 表层 回 波 丰 ，(C) 考虑 到 来 自 星 下 点 外 杂 波 的 次 表层 探测 雷达 


原理 示意 图 1。 


Fig.1 (A) Operation principle of the space-based subsurface investigation radar”! ;(B) Illustration of distribution of 


surface and subsurface echo in time/range. Strong surface echo masking a weaker subsurface return; 


(C)Operation principle of the space-based subsurface investigation radar in view of off-nadir surface scatters 03]. 


由 于 雷达 安装 在 环绕 器 上 ， 与 探测 目标 的 距离 至 少 为 轨道 高 度 ， 按 照 雷 达 分 辨 理论 ， 
雷达 的 方位 向 分 辩 率 为 H .A (2D)， 在 次 表层 探测 雷达 的 运行 的 情境 下 方位 向 分 辨 率 将 达到 
BOP TAK, 要 提高 雷达 的 方位 向 分 辨 率 只 能 通过 增 大 雷达 天 线 的 口径 , 而 在 太空 中 这 种 要 求 
是 不 切实 际 的 ， 因 此 次 表层 探测 雷达 都 采用 了 合成 孔径 技术 ， 将 在 不 同时 刻 发 射 、 接 收 的 电 
伺 波 等 效 为 在 同一 时 刻 发 射 、 接 收 ， 即 合成 为 一 个 大 孔径 天 线 进行 发 射 和 接收 。 在 这 种 技术 
的 帮助 下 ,合成 孔径 大 小 可 以 视 为 V24， 将 其 代入 方位 向 分 辩 率 公式 ,将 大 大 提高 方位 向 分 
辨 率 ， 最 优 方 位 向 分 辨 率 可 以 达到 VX 互 /(2V2)" ,因此 轨道 器 次 表层 探测 雷达 在 数据 处 理 
中 必须 进行 方位 向 处 理 来 实现 合成 孔径 。 

电磁 波 从 发 射 机 通过 天 线 发 射 ， 穿 越 电 离 层 ， 被 目标 反射 再 次 穿越 电离 层 ， 最 终 被 雷 
达 接 收 ， 这 个 过 程 在 仿真 中 可 以 认为 是 理想 的 。 实 际 上 ， 天 线 和 接收 通道 会 对 电磁 波幅 度 和 
相位 造成 影响 , 电离 层 造成 的 相位 失真 也 严重 影响 了 成 像 质 量 , 为 了 保证 雷达 的 成 像 精 度 和 


一 


ninaXiv 合 作 期 刊 


142 ”次 表层 探测 的 效果 ， 必 须 对 上 述 过 程 造 成 的 幅 相 失 真 进行 校正 。 


143 
144 1.3 轨道 器 次 表层 探测 雷达 数据 处 理 流程 与 方法 
145 合成 孔径 雷达 数据 处 理 过 程 一 般 可 以 划分 为 距离 向 处 理 和 方位 向 处 理 (或 多 普 勒 处 理 ) ， 


146 ”在 MARSIS 与 SHARAD 的 数据 处 理 流程 中 ，MARSIS 较 常用 的 探测 模式 SS3 (subsurface 3, 
147 ”次 表层 探测 3 模式 ) 选择 在 星 上 完成 方位 向 处 理 后 下 传 ， 再 在 地 面 完 成 距离 向 处 理 ; 
148 ”SHARAD 的 原始 数据 在 星 上 完成 相关 预 全 加 后 ， 先 进行 距离 向 处 理 ， 再 进行 方位 向 处 理 。 
149 ”其 中 还 包括 对 电离 层 带 来 影响 的 校正 ， 使 用 的 方法 也 有 不 同 。 

150 ”1.3.1 距离 向 处 理 

151 合成 孔径 雷达 数据 处 理 中 ， 距 离 向 处 理 的 操作 都 为 脉冲 压缩 ， 将 接收 到 的 回 波 信号 通 
152 ”过 匹配 滤波 器 , 将 回 波 信号 与 发 射 信号 进行 相关 处 理 , 这 种 操作 能 将 时 间 上 分 散 的 信号 能 量 
153 ”集中 到 一 个 时 间 点 附近 , 使 得 来 自 地 面 或 次 表层 界面 的 回 波 时 间 宽 度 变 窄 ,能量 也 被 压缩 到 
154 ”该 时 间 点 。MARSIS 与 SHARAD 进 行 脉冲 压缩 的 操作 在 细节 上 有 一 些 不 同 ，MARSIS 的 脉冲 
155 ”压缩 流程 为 : 


p 
V | 


0 


=p 
/ 
g 


= 156 1) 对 回 波 信 号 进行 一 次 快速 傅 里 叶 变 换 CFFTO ， 点 数 为 512; 
157 2) 对 作为 参考 信号 的 线性 调频 信号 进行 傅 里 叶 变 换 , 求 得 其 频谱 , 再 对 其 进行 共 轿 变换 ; 
158 3) 将 1) 、2) 中 得 到 的 回 波 信 号 频谱 与 参考 信号 频谱 共生 相 乘 ， 进 行 反 传 里 叶 变换 得 到 


159 PEAR", 
于 天 线 与 收发 通道 、 电 离 层 造成 的 幅 相 变化 ， 如 果 使 用 理想 的 线性 调频 信号 作为 参 
161 ” 考 信号 ， 则 其 与 回 波 信和 号 的 相位 无 法 完全 匹配 ， 将 影响 脉冲 压缩 的 效果 。 因 此 MARSIS 使 用 
162 ”的 参考 信号 是 由 理想 线性 调频 信号 与 在 仪器 定 标 中 测 得 的 幅 相 失真 推导 得 到 的 , 并 根据 电离 
163 ” 层 失 真 对 参考 信号 进一步 调整 。 此 外 ， 为 了 降低 脉 压 后 出 现 的 旁 兴 ,MARSIS 使 用 了 汉 宁 窗 
164 《hanning〉 对 参考 信号 进行 加 窗 处 理 ， 加 窗 可 以 降低 脉 压 结果 的 旁 渡 ,但 也 会 导致 主流 展 
165 
166 SHARAD 探 测 数据 下 传 后 通过 两 种 不 同 的 数据 处 理 流 程 完 成 距离 向 脉冲 压缩 ， 第 一 种 
167 ”流程 的 参考 信号 基于 天 线 匹 配 网 络 定 标 时 测 得 的 信号 幅度 , 测 得 的 定 标 数据 为 接收 机 与 发 射 
168 ”机 温度 的 函数 。 使 用 定 标 数 据 得 到 的 参考 信号 完成 脉 压 后 ， 使 用 相位 梯度 自 聚 焦 算法 
169 (phase-gradient autofocusing, PGA) 进行 电离 层 校正 。 而 第 二 种 流程 选择 使 用 线性 调频 信 
170 ”号 频率 分 量 的 等 幅 模 型 , 结果 是 脉 压 后 旁 准 将 出 现 一 种 不 对 称 偏 置 , 时 延 较 大 的 旁 汶 幅度 较 
171 K, 选择 这 种 方法 是 因为 其 可 以 在 脉 压 完 成 后 保持 两 倍 的 过 采样 率 , 有 利于 在 接 下 来 的 处 理 
172 ”中 进行 插值 。 两 种 处 理 流程 都 选择 了 汉 宁 窗 对 数据 进行 加 窗 以 减 小 旁 久 。 在 进行 电离 层 校 正 
173 ”方面 , 第 二 种 流程 使 用 了 一 种 频率 决定 的 相位 误差 模型 , 这 种 模型 给 出 的 校正 项 与 总 电子 含 
174  i& (TEC) 接近 线性 相关 ， 因 此 与 来 自 地 表 和 次 表层 回 波 信 号 时 延 的 改变 量 相关 ， 这 种 技术 
175 ”适用 于 探测 时 太阳 天 顶 角 小 于 100” 获 得 的 图 像 。 
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1. 3.2 方位 向 处 理 

方位 向 处 理 是 合成 孔径 雷达 技术 的 关键 与 难点 所 在 ， 不 同 的 处 理 方法 需要 的 处 理 量 
不 同 。 轨 道 器 次 表层 探测 雷达 的 方位 向 处 理 方法 与 普通 合成 孔径 雷达 的 方位 向 处 理 方法 在 原 
理 上 是 相同 的 , 因此 轨道 器 次 表层 探测 雷达 数据 处 理 使 用 的 方位 向 处 理 算法 都 是 在 常见 的 合 
成 孔径 雷达 算法 的 基础 上 产生 的 。 比 较 常 见 的 合成 孔径 雷达 方位 向 处 理 算 法 包括 CS 算法 
(chirp scaling algorithm)0s、ok 算 法 上 由]、RD 算 法 Crange-doppler algorithm ,距离 -多 普 勒 算法 ) 
99、BP 算 法 (back-projection algorithm, 后 向 投影 算法 ) 四 等 ， 如 果 对 数据 进行 非 聚 焦 成 像 ， 
则 可 以 选择 多 普 勒 波束 锐 化 算法 中" 等 。 

RD 算法 在 上 世界 70 年 代 末 就 已 经 被 用 于 合成 孔径 雷达 成 像 ， 至 今 仍 在 广泛 使 用 。 RD 
法 的 基本 实现 方法 为 对 雷达 回 波 数据 矩 阵 进行 距离 向 脉冲 压缩 ,然后 对 数据 进行 方位 向 傅 里 
叶 变 换 ， 将 数据 转换 到 方位 向 频 域 ， 并 用 sinc 函 数 进行 距离 插值 完成 距离 徙 动 校正 ， 将 距离 
徙 动 曲线 拉 直 到 与 方位 向 频率 轴 相 平行 的 方向 ,最 后 进行 方位 向 压缩 完成 成 像 。 这 种 算法 通 
过 距离 向 与 方位 向 的 频 域 处 理 ， 所 有 操作 都 在 一 维 频 域 内 进行 ， 简 便 高 效 ， 但 由 于 校正 距离 
徙 动 所 需 的 插值 操作 ， 其 对 长 孔径 SAR 的 处 理 精 度 有 限 ， 且 运算 量 较 大 。CS 算 法 避免 了 插 
值 操作 ， 基 于 chirp scaling 原 理 ， 通 过 相位 相 乘 代替 时 域 插值 以 完成 距离 徙 动 校正 ， 并 将 所 
有 信和 号 的 距离 徙 动 分 为 补 余 距 离 徙 动 与 一 致 距离 徙 动 ， 即 将 距离 徙 动 分 解 出 一 致 的 一 部 分 
(BUSES) ， 二 者 的 差 为 补 余 距 离 徙 动 。CS 算 法 的 处 理 流 程 首先 将 数据 进行 方位 向 
傅 里 叶 变换 ,进行 相位 相 乘 补偿 补 余 距 离 徙 动 , 然后 进行 距离 向 傅 里 叶 变换 将 数据 变换 到 二 
维 频 域 ， 与 参考 函数 进行 相位 相 乘 ， 完 成 距离 压缩 、 二 次 距离 压缩 以 及 一 致 距离 徙 动 校正 ， 
接 下 来 进行 距离 向 傅 里 叶 道 变换、 方位 向 压缩 、 方 位 向 傅 里 叶 逆 变 换 ， 完 成 成 像 。 还 有 一 种 
ok 算法 在 处 理 宽 孔 径 或 大 斜视 角 数 据 上 具有 独特 的 优势 ， 这 种 算法 使 用 二 维 傅 里 叶 变 换 将 
数据 变换 到 二 维 频 域 ， 并 乘 上 参考 函数 ， 使 参考 距离 上 的 目标 得 到 全 聚焦 ， 其 他 目标 得 到 部 
分 聚焦 , 接 下 来 进行 Stolt 插 值 完成 其 他 目标 的 聚焦 ,将 数据 变换 回 时 域 完成 处 理 中 。 这些 算 
法 有 各 自 的 优势 与 缺陷 , 适应 于 不 同 的 雷达 成 像 场景 , 因此 算法 的 选择 需要 结合 雷达 数据 的 
具体 情况 进行 分 析 。 

SHARAD 选 择 了 CS 算法 ， 将 方位 向 上 回 波 的 多 普 勒 频 移 进行 分 解 ， 将 距离 徙 动 划分 为 
不 随 距 离 变 化 的 一 致 距离 徙 动 与 随 距 离 变 化 的 补 余 距离 徙 动 ， 在 处 理 过 程 中 分 别 进行 校正 。 
这 种 算法 只 需要 快速 傅 里 叶 变 换 与 复 乘 就 能 完成 中 ， 计 算 量 较 小 ， 使 用 简单 。 具 体 的 处 理 
流程 为 : 1) 计算 距离 徙 动 ， 将 每 道 数 据 的 每 个 时 延 位 置 点 相对 物理 孔径 中 心 的 距离 进行 排 
列 ; 20 方位 压缩 ， 用 环绕 器 的 轨道 高 度 ， 径 向 速度 ， 切 向 速度 确定 每 个 回 波 数据 的 差分 相 
位 ， 校 正 相 位 可 以 表示 为 整个 带宽 的 频率 的 函数 ， 但 实际 操作 中 使 用 最 大 的 工作 频率 (25 
MHz) 可 以 达到 最 佳 的 聚焦 效果 ; 3) 完成 相位 校正 后 ,对 完成 距离 徙 动 校正 的 “等 时 延 线 ” 
进行 傅 里 叶 变 换 ， 取 所 得 频谱 的 幅 值 即 为 一 列 雷达 图 像 ， 使 用 汉 宁 窗 在 方位 向 上 压低 旁 淤 ， 
在 顺 轨 方向 约 每 460 m. 〈1 个 度 间 隔 ) 上 有 128 列 图 像 。 如 果 多 普 勒 频 域 的 单位 频 宽 比 一 个 频 
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率 分 辨 单元 更 宽 ， 则 可 以 进行 多 视 处 理 。 

MARSIS 选 择 的 合成 孔径 处 理 算 法 与 SHARAD 的 CS 算法 有 巨大 区 别 ， 这 种 算法 被 称 为 
多 普 勒 波束 锐 化 , 在 完成 对 轨道 器 垂直 方向 的 运动 补偿 后 进行 ,可 以 补偿 轨道 器 除 垂 直方 向 
外 其 他 方向 运动 带 来 的 相位 。 这 是 一 种 非 聚 焦 SAR 技 术 ， 获 得 的 图 像 方位 分 辩 率 较 低 ， 但 计 
算 量 相对 较 小 。 人 处 理 方 法 为 将 一 帧 (frame〉 内 的 所 有 同 频率 的 回 波 做 又 加 处 理 ， 对 其 相位 
进行 调整 完成 相关 性 加 。 一 “ 帧 ” (frame〉 中 包含 了 在 一 个 合成 孔径 时 间接 收 到 的 回 波 
信和 号， 由 工作 频率 /波长 以 及 雷达 高 度 决定 。 经 过 欠 加 的 回 波 信 号 为 一 个 多 普 勒 滤波 器 ， 多 
普 勒 滤波 器 可 以 在 接 下 来 的 处 理 中 由 地 面 控 制 选择 前 表面 反射 法 或 对 比 法 , 对 孔径 中 心 的 多 
普 勒 滤波 器 第 0 多 普 勒 滤波 器 进行 电离 层 相 位 校正 ， 校 正 的 结果 将 在 用 于 更 正 距离 压缩 
中 的 参考 函数 。 


full Rate Processing ae Rate Processing 


图 2 MARSIS 多 普 勒 (方位 向 处 理 与 距离 向 处 理 流程 图 站。 


Fig.2 MARSIS range processing and doppler processing flow diagram !?l. 


而 我 国 首次 火星 探测 任务 “天 问 1 号 ”的 轨道 器 次 表层 探测 雷达 (MOSIR) 在 数据 处 理 
流程 设计 的 过 程 中 , 为 实现 星 上 实时 成 像 处 理 , 初步 使 用 了 BP 算法 作为 方位 问 处 理 的 算法 ， 
而 对 于 下 传 的 原始 数据 ,在 地 面 采 用 何 种 方法 进行 处 理 ， 则 需要 结合 数据 的 实际 情况 ,选择 
成 像 效 果 最 好 、 距 离 徙 动 计算 最 准确 的 方法 。 

1.3.3 电离 层 校正 

MARSIS 次 表层 探测 的 工作 模式 与 火星 电离 层 等 离子 体 频 率 非 常 接近 09, 因此 电离 层 对 
回 波 时 延 以 及 幅度 、 相 位 的 影响 是 不 可 忽略 的 。 在 最 初 的 处 理 中 , 研究 人 员 提 出 了 两 种 校正 
算法 , 分 别 为 前 表面 反射 法 与 对 比 法 , 文 [17] 等 文献 中 更 多 地 使 用 对 比 法 来 完成 电离 层 校正 ， 
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校正 效果 见 图 3。 对 比 法 通过 对 电离 层 相 位 误差 的 泰勒 展开 式 中 影响 较 大 的 系数 进行 估计 以 


完成 校正 ， 并 求 得 电离 层 电 子 总 含量 (TEC) 


， 对 回流 进行 迭代 脉 压 ， 在 每 次 脉 压 时 对 回 波 


这 号 施加 估计 的 相位 误差 校正 项 , 估计 相位 误差 使 用 的 参数 步 长 由 发 射 的 线性 调频 信号 以 及 
回 波 信号 的 调频 率 决 定 ， 直 到 校正 精度 满足 要 求 ， 校正 后 的 回 波 信号 能 量 集 中 程度 最 高 ， 此 
时 根据 完成 校正 所 得 的 泰勒 展开 式 系数 可 以 计算 出 TBCP。 


图 3 未 校正 电离 层 失真 的 MARSIS 第 5017 轨 雷 达 图 


f& (上 图 ) 与 使 用 对 比 法 进行 电离 层 校正 后 (下 图 ) 


fpe m. 


Fig.3 Radargram for orbit 5017 before(top) and after(bottom) use of the contrast method to correct ionospheric 


distortion! 


17]. 


SHARAD 的 工作 频段 较 高 ， 而 火星 电离 层 的 等 离子 体 频 率 一 般 为 几 个 MHz 数 量 ， 频 率 
高 于 电离 层 等 离子 体 频 率 的 电磁 波 在 穿越 电磁 波 时 都 会 发 生 衰减 , 信号 传播 时 间 将 出 现 群 延 


迟 。 在 SHARAD 的 星 上 处 理 过 程 中 ， 相 位 梯度 自 聚 焦 算 法 PGA) 用 于 进行 校正 电离 层 带 


来 的 相位 误差 ， 在 MARSIS 的 后 期 数据 处 理 中 


， 这 种 算法 也 被 作为 对 比 法 等 的 补充 中。 这 种 


算法 可 以 同时 对 电离 层 以 及 其 他 来 源 造 成 的 相位 误差 进行 估计 ， 上 其 体操 作 需 要 选取 特 显 点 ， 
即 含有 较 高 能 量 且 未 经 补偿 的 帧 ,利用 散 焦 图 像 中 相位 误差 的 元 余 信 息 对 相位 误差 的 梯度 进 
行 鲁 棒 估计 ， 一般 需要 经 过 圆 移 ( 将 特 显 点 移 至 方位 中 心 ) ， 加 窗 ， 相 位 梯度 估计 ， 最 后 进 
行 近 代 校 正 ,完成 整个 校正 过 程 ,需要 注意 的 是 , PGA 校正 的 效果 与 特 显 点 选择 有 很 大 关系 ， 

如 果 选 取 的 特 显 点 帧 中 , 相位 误差 泰勒 展开 式 中 影响 较 大 的 系数 不 稳定 , 则 PGA 的 校正 效果 


将 被 削弱 。 


1.4 MARSIS & SHARAD 成 像 结 果 对 比 
从 MARSIS 与 SHARAD 的 工作 参数 以 及 型 


LE 论 探测 能 力 中 ， 可 以 看 出 MARSIS 的 工作 频率 


低 ， 带 宽窄 ， 波 长 长 ， 且 合成 孔径 处 理 采 用 的 是 非 聚 焦 的 多 普 勒 波束 锐 化 算法 ，SHARAD 
的 工作 频率 较 高 , 带宽 宽 , 波长 短 , 使 用 chirp scaling 算 法 进行 合成 孔径 处 理 。 因 此 , MARSIS 


的 探测 深度 较 深 , 可 以 探测 到 距 地 表 3km 的 次 


表层 界面 SHARAD 接收 到 的 来 自 相同 深度 的 


回 波 太 弱 以 至 于 被 背景 噪声 淹没 无 法 探测 ,但 MARSIS 的 探测 精度 无 论 是 距离 分 状 率 还 是 顺 
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轨 分 辩 率 还 是 交 轨 分 辨 率 都 明显 弱 于 SHARAD， 尤 其 是 距离 分 辨 率 ，SHARAD 的 理论 距离 
分 辨 率 只 有 MARSIS 的 十 分 之 一 。MARSIS 与 SHARAD 在 相同 区 域 的 探测 图 像 可 以 作为 互补 
资料 ， 根 据 MARSIS 探 测 数据 在 较 大 尺度 上 绘制 探测 区 域 的 图 像 ， 确 定 有 精细 探测 价值 的 目 
标 区 域 ， 再 利用 SHARAD 的 数据 对 该 区 域 火星 地 表 以 下 数 百 米 可 能 存在 的 次 表层 结构 进行 
T A HERE RI o 

文 [15] 对 相近 区 域 的 MARSIS 与 SHARAD 成 像 进行 了 比较 ，MARSIS 28 3749 9), 55 
SHARAD 第 2026 轨 探测 区 域 非常 接近 ， 可 以 比较 两 部 雷达 对 火星 相同 区 域 成 像 的 效果 ， 如 
图 4 所 示 吕 。 可 以 看 到 SHARAD 对 距 地 表 回 波 约 7 us ( 约 数 百 米 深 处 ) 次 表层 界面 成 像 清晰 ， 
能 反映 出 次 表层 较 精细 的 结构 ， 但 探测 深度 较 MARSIS 较 浅 ， 实 际 成 像 质 量 与 理论 计算 的 探 
测 精度 相符 。 


SHARAD RADAR GRAM: Detail 


图 4.MARSIS 第 3749 轨 与 SHARAD 第 2026 轨 探 测 区 域 重合 部 分 图 像 比较 (9。 


Fig.4 Comparison between MARSIS radargram orbit 3749 and SHARAD radargram orbit 2026 on the same Mars 


area ll. 


1.5 研究 现状 

MARSIS 与 SHARAD 运 行 近 20 年 来 , 基本 实现 了 对 火星 全 球 的 覆盖 ,截至 2019 年 仍 在 发 
布 数据 , 研究 人 员 根 据 其 发 布 的 数据 , 完成 数据 成 像 处 理 后 , 使 用 搭载 在 “火星 全 球 勘探 者 ” 
号 探测 器 上 的 火星 轨道 激光 高 度 计 (MOLA) 获得 的 高 精度 火星 表面 高 程 数据 (空间 分 辨 率 
接近 1”， 绝 对 精度 达到 13m! D ， 建 立 探测 区 域 火 星 表面 地 形 模 型 ， 仿 真 雷达 接收 到 的 表 
面 杂 波 ， 短 除 与 次 表层 回 波 发 生 混 杂 的 表面 杂 波 中 I。 获得 次 表层 回 波 后， 结合 周边 的 地 
形 情 况 ， 设 计 多 种 方法 反 演 出 火星 表层 、 次 表层 物质 的 介 电 常 数 、 衰 减 系 数 、 次 表层 界面 深 
度 等 科学 信息 ， 推 测 火星 表层 、 次 表层 物质 成 分 ， 获 得 了 众多 高 价值 的 研究 成 果 ， 以 此 为 依 
据 提 出 了 关于 火星 的 气候 、 地 质 、 轨 道 等 历史 演化 的 一 系列 推论 。 

火星 北极 高 原 是 MARSIS 的 首 个 探测 区 域 , 覆盖 了 火星 北极 周边 直径 约 600 公 里 的 面积 ， 
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主要 成 分 为 水 冰 与 少量 尘埃 ， 由 三 个 地 质 单元 组 成 :北方 残余 冰 盖 《NRIC〉、 北 极 层 状 沉 
$R (NPLD) 、 基 底层 (the basal unit) 。 文 [26] 利 用 MARSIS 早 期 中 心 频率 为 3 MHz 和 5 MHz 
在 第 1855 轨 所 得 探测 数据 , 获得 了 从 北部 平原 穿 过 北极 层 状 沉积 边缘 进入 北极 层 状 沉积 区 域 
的 雷达 图 像 , 雷达 图 像 显 示 表 面 反 射 回 波 分 裂 成 了 两 个 强 反 射 回 波 , 使 用 MOLA 数 据 对 探测 
区 域 表 面 回 波 进行 仿真 , 排除 了 较 低 的 反射 回 波 是 由 表面 杂 波 导致 的 可 能 性 , 即 说 明 该 反射 
回 波 来 自 次 表层 分 层 界面 ,两 个 反射 回 波 之 间 的 时 延 以 及 相对 回 波 强度 说 明基 底 上 禾 物 质 的 
介 电 常 数 与 损耗 角 正 切 与 纯 水 冰 的 介 电 特性 相符 , 用 纯 水 冰 的 介 电 常 数 (~3) 将 时 延 转换 为 
深度 ， 发 现 反 射 体 的 深度 与 北极 层 状 沉积 高 于 周围 地 区 的 高 度 大 致 相等 ， 为 1.8km， 说 明 北 
极 层 状 沉积 下 的 次 表层 界面 仍 与 周边 地 区 的 平原 平 齐 , 没有 向 下 弯曲 , 证 明了 非常 厚 的 弹性 
岩石 圈 的 存在 ， 还 表明 壳 幅 温度 梯度 较 低 。 


a 


g 
Elevation, m o 


-5000 -4500 -4000 -3500  -3000 -2500 


图 5 CAO 第 1855 轨 MARSIS 探测 数据 (轨道 从 周边 地 区 进入 NPLD〉; (B) 使 用 MOLA 数据 模拟 的 
MARSIS 探测 到 的 表面 回 波 (包括 星 下 点 回 波 以 及 星 下 点 外 的 杂 波 ) ; CC) 沿 火 星 快 车 (Mars Express? 
探测 器 星 下 点 轨迹 绘制 的 MOLA 数据 地 形 图 2 。 


Fig.5 (A) MARSIS data from orbit 1855 as it crossed the margin of the NPLD and came into the NPLD.(B) 
Simulated MARSIS data if echoes are only from the surface(nadir and off-nadir clutter). (C) MOLA topography 
along the ground track(red line) of the Mars Express sub-satellite point), 

南极 高 原 与 北极 高 原 相 似 ， 也 由 三 个 主要 地 质 单元 组 成 : 南方 残余 冰 羡 CSRIC) 、 南 
极 层 状 沉积 (SPLD)、 阿 詹 泰 山 浓 构造 (DAF) ,南方 残 余 冰 羡 表 面 履 盖 着 数 米 厚 的 干冰 ， 
具有 较 高 的 反照 率 ， 南 极 层 状 沉积 的 结构 与 北极 层 状 沉积 相似 。 

文 [27] 根 据 MARSIS 第 2753 轨 的 探测 数据 ， 发 现 星 下 点 在 通过 南极 层 状 沉积 边缘 时 ， 也 
出 现 了 表面 回 波 分 裂 成 两 个 反射 回 波 的 现象 , 且 较 低 的 回 波 使 用 水 冰 的 介 电 常数 将 时 延 换算 
为 深度 后 , 次 表层 分 层 界面 符合 南极 层 状 沉积 周边 地 区 表面 地 形 的 延伸 , 该 界面 被 认为 是 南 
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极 层 状 沉积 的 富 冰 材料 与 以 岩 眉 为 主要 成 分 的 基底 的 分 界 ,MARSIS 在 南极 层 状 沉积 的 大 部 
分 区 域 都 探测 到 了 这 个 界面 。 从 图 6 中 可 以 看 出 次 表层 界面 回 波 的 幅度 很 强 ， 甚 至 在 某 些 区 
域 (如 第 2682 轨 ) 次 表层 界面 回 波 强 于 表面 回 波 ,这 说 明 南 极 层 状 沉积 的 成 分 具有 极 低 的 衰 
减 系数 ， 损 耗 角 正切 约 为 0.001~0.005， 可 以 由 此 推测 出 南极 层 状 沉积 中 水 冰 的 含 侍 量 在 0 到 
10% 之 间 。 
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图 6 (A) MARSIS 第 2753 轨 探 测 数据 ; CB) 第 2753 轨 星 下 点 轨迹 周围 MOLA 数据 地 形 图 ;(C) MARSIS 
第 2682 PIRMA; D) 第 2682 轨 星 下 点 轨迹 周围 MOLA 数据 地 形 图 ， Œ) 第 2682 轨 探 测 数 据 ， 
MOLA 测 得 表面 轮 廊 (黑色 实 线 与 MARSIS 测 得 的 基底 形状 〈 蓝 色 ) 971, 


Fig.6 (A) MARSIS data from orbit 2753. (B) MOLA topography along the ground track of orbit 2753. (C) 


MARSIS data from orbit 2682. (D) MOLA topography along the ground track of orbit 2682. (E) MARSIS data 
from orbit 2682. MOLA surface elevations (black line) and MARSIS measured basal elevations (blue symbols) 中 
ZEITMARSIS XY FA B DURAS DU ARR, SC[277 43s HI nii I PRR 
MARSIS 探 测 数据 与 相应 的 MOLA 高 分 辨 率 火 星 表 面 地 形 高 程 数 据 ，MOLA 表 面 高 程 数据 减 
去 经 过 插值 的 MARSIS 测 得 的 基底 高 程 数 据 ， 得 到 对 应 位 置 的 南极 层 状 沉积 的 厚度 ， 并 绘制 
厚度 分 布 图 〈 见 图 7) 。 
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图 7 基于 MARSIS 测 得 数 ] 


与 MOLA 表面 高 程 数据 绘制 火星 南极 层 状 沉积 厚度 分 布 图 ?7。 


Fig.7 Map of the SPLD thickness, based on MARSIS measurements and MOLA surface topography °”. 


MARSIS i FR fg XE CI I BLA DE A PARRA,  XC[28] T MARSIS 1E FA BR r FE 
DC PA193* E, 81^ S$ 为 中 心 ， 宽 度 为 200km， 大 部 分 地 区 
回 波 十 分 强烈 , 甚至 强 
绘制 成 归 一 化 后 


地 区 进行 探测 获得 的 数据 ， 探 测 
回 波 都 较 弱 , (AA 


的 基底 反射 


3 


将 次 表层 回 波 的 强度 归 一 化 到 表面 回 波 的 强度 , } 
Kl, 分布 图 显示 ， 在 193”E，81”S 附 近 存 在 一 个 20km 宽 的 次 表层 回 波 异 常 区 
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下 可 能 存在 液态 水 ， 而 水 中 含有 鳞 、 钙 、 钠 等 的 高 氧 酸 盐 


原因 。 


文 [29] 使 月 


进行 了 研究 ， 


数 
体积 的 估计 范 


电 常 数 实 部 为 2.9 土 0.4, E 


结果 暗示 梅 杜 莎 槽 沟 沉 积 


E 某 些 区 域 , 基底 反射 
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组 归 一 化 后 基底 回 波 能 量 


与 基底 介 电 常 数 关 系 


的 可 能 


HH 线 的 包 络 


性 ， 


分 析 仿 真 结果 得 到 异常 


区 域 


的 介 电 常数 取 值 大 于 15, 这 一 结论 拟 


出 了 


月 MARSIS 对 火星 经 度 130"” ~240° E 赤 道 地 
梅 杜 水 槽 沟 构 造 可 能 是 火星 上 最 年 轻 的 地 表 沉 积 
105km2。MARSIS 探 测 到 了 梅 杜 莎 模 沟 沉积 物 与 下 方 地 形 构造 的 次 表层 界面 


据 揭 示 了 梅 杜 莎 槽 沟 构 造 沉积 物 的 厚度 与 介 电 特 性 
Fl (1.4~1.9X10°km?) 。 由 探测 数据 


合理 解释 : 该 异常 


口 
是 其 在 低温 环境 下 仍 能 保持 液态 的 


区 域 的 层 状 沉积 


区 的 探测 数据 ， 对 梅 杜 莎 槽 沟 构 造 
区 域 之 一 ， 面 积 约 为 2.1 X 
(图 10) ， 探 测 


E, FABRA SEF ATER 


PU BOISE EZ H 


物 


ua 


物 ， 如 果 沉 积 物 为 后 者 ， 则 其 中 
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日 


沟 沉积 物 


bs 梅 杜 莎 权 沟 沉积 物 的 介 
次 表层 回 波 的 损耗 计算 出 衰减 系数 为 ~0.0048 土 0.0024dB/m。 上 述 
能 为 干燥 的 疏松 物质 或 水 冰 与 高 介 电 常数 的 其 他 物质 的 混合 


的 含 侍 量 将 高 于 南极 层 状 沉 积 的 10%。 
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图 8 MARSIS 第 10737 轨 在 南极 高 原 获得 的 探测 数据 ， 可 见 部 分 地 区 基 岩 回 波 极其 强烈 ， 
形成 高 亮 区 域 p%。 


Fig.8 MARSIS data from orbit 10737 on the Planum Australe. Strong basal reflections can be seen at some 


locations, where show as highlight area 8), 
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图 9 〈 左 ) MARSIS 中 心 频 率 4MHz 下 对 异常 区 域 地 质 分 层 模 型 进行 仿真 的 结果 ， 蓝 色 阴 影 区 域 为 归 一 化 
后 基底 回 波 能 量 与 基底 介 电 常 数 关 系 的 


| 线 包 络 ， 黑 色 直 线 为 实际 探测 数据 中 基底 回 波 能 量 中 值 ，( 右 )》 
或 内 外 基底 介 电 常数 取 值 分 布 ?™。 


uu 


异常 区 


Fig.9 (Left) Results of the electromagnetic simulation of layers model of the bright reflection area computed at 4 
MHz The blue shaded area is the envelope of all curves between normalized basal echo power and basal 
permittivity, and the black horizontal line is the median normalized basal echo power at 4 MHz from the 

observations; (Right) Basal permittivity distributions inside (blue) and outside (brown) 


the bright reflection area '*), 
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图 10 (A) MARSIS 第 2896 轨 探 测 数 据 ; 
测 数据 ; (D)MARSIS 第 3824 轨 探 测 数 ] 


据 ; 


(B) MARSIS 第 


E; (EMARSIS 第 4117 gud Xi] 


4011 轨 探 测 
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以 上 探测 数据 均 在 梅 杜 莎 槽 沟 构 造 区 域 来 集 。 


数据 ，(C)MARSIS 第 3868 轨 探 


fi; (F)MARSIS 第 3996 轨 探 测 数 


Fig.10 (A)MARSIS data from orbit 2896; (B)MARSIS data from orbit 4011; (C)MARSIS data from orbit 


3868;(D)MARSIS data from orbit 3824;(E)MARSIS data from orbit 4117;(F)MARSIS data from orbit 3996; these 


MARSIS orbit tracks were measured in the Medusae Fossae Formation 


SHARAD 由 于 其 较 高 的 频率 和 较 宽 


[29] 


穿 透 深度 较 浅 ， 但 深度 分 辩 率 较 高 ， 其 对 


南极 高 原 、 北 极 高 原 以 及 其 他 地 区 可 以 进行 更 为 精细 的 探测 ， 获 取 了 一 些 MARSIS 无 法 探测 


到 的 


原 的 
为 两 


信息 。 


文 [30] 利 用 SHARAD 在 火星 


a 


个 时 延 不 同 的 回 波 〈 图 11) ， 


Esser ep RED 
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p mu m Lu RS YP ARE 
超过 100 轨 探测 数据 ， 发 现 星 下 点 在 穿 过 河谷 平原 经 过 舌 状 岩 屑 坡 时 出 
HMOLA Zi € Bi; Ex £5 RN 


现 表 面 加 


波 分 列 


KM 


ARDS BL, VA BC ORAS JE IRAE. EEA TARA TP E A FL CS EA) 


为 3， 


都 有 分 布 。 文 [31] 使 月 
测 数据 ， 绘 制 了 多 幅 舌 状 岩 居 
对 比 以 及 多 轨 探 测 数据 间 的 对 比 , 


表明 


越 介质 造成 的 衰减 ， 以 及 表面 回 波 与 次 表层 回 
水 冰 相 符 ， 且 证 明 血 状 岩 眉 坡 下 掩埋 的 冰川 含 尘 量 
测 得 雷达 波 速 与 水 冰 中 的 雷达 波 速 类 似 ， 与 MARSIS 在 极地 层 状 沉积 的 观测 结果 相符 ， 


与 纯 


及 其 


次 表层 介质 十 分 均一 ,上 且 


推测 得 衰减 率 为 ~2dB/h sec， 与 较 纯净 
这 种 舌 状 岩 届 坡 (LDAs) 不 仅 出 现在 北半球 中 纬度 地 
SHARAD 在 火星 南半球 中 纬度 的 海 勒 
坡 的 次 表层 雷达 图 
上 定 了 表面 以 下 的 反射 
.大 型 岩石 碎片 不 是 其 中 的 主要 组 分 。 洲 
波 的 关联 ， 测 得 衰减 系数 最 大 为 ~10dB/km， 
不 可 能 高 于 10%。 较 弱 的 衰减 率 ， 以 


象 ( 图 12 


这 一 


结果 指出 舌 状 岩 眉 坡 内 部 的 3 


要 成 分 为 水 冰 , 印记 


的 水 冰 或 水 冰 与 少量 其 他 物质 的 混合 物 相 符 。 
区 ， 在 火星 南北 半球 30”-~60” 
斯 撞击 坑 东 部 边缘 附近 的 多 组 探 


) ， 经 过 与 MOLA 模 拟 仿真 数据 的 
回 波 为 次 表层 界面 回 波 ， 分析 图 像 
| 量 探测 数据 中 雷达 波 穿 
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Fig.11 SHARAD radargram along the ground track of lobate debris aprons 
星 北 半球 中 纬度 的 乌托邦 低 原 直径 约 为 3300km， 多 见 多 
究 人 员 根 据 对 乌托邦 低 原 进行 的 形 貌 分 析 猜 测 该 区 域 近 地 表 含 有 


zn 
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图 11 SHARAD Æ PAFITA SIN RS DUIS 。 


[30] 


边 形 与 扇形 洼地 


MOLA 测 得 的 火星 表面 数字 高 程 模型 ， 利 
常数 为 2.8 土 0.8， 符 合 水 冰 、 石 质 物质 以 及 
E50%-85% 之 间 ， 石 质 材料 低 于 30%， 了 筷 隙 空洞 则 占 
到 了 体积 的 15%-50%， 而 水 冰 体 积 达到 了 8400-14300km， 乌 托 邦 低 原 存在 水 冰 的 说 法 得 到 


SHARAD 探 测 数 据 
体 深度 ， 拟 合 后 估计 近 次 表层 物质 的 介 昌 
EE 特性， 估计 其 中 水 冰 含 量 


J PEW T RR 


大 量 水 冰 。 文 [32] 运 用 


雷达 探测 数据 、 火 星 气候 模型 以 及 表 再 
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图 12 (A) MOLA 数据 仿真 得 到 的 表面 杂 波 ; 


形 貌 分 析 的 支持 。 
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E 实 存在 的 回 波 ， CC) 将 原始 探测 数据 i 


加 像 ， 假 设 舌 状 岩层 坡 表面 下 的 成 分 为 水 冰 5 1。 
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Fig.12 (A)Simulated surface clutters with MOLA data. (B)SHARAD data from orbit 6830, Vertical arrows identify 


echoes not consistent with surface clutter simulation and also confirmed in adjacent tracks. (C)Radar data 


(31] 


converted into depth assuming a water-ice composition ^. 


SHARAD 在 南极 高 原 、 北 极 高 
候 变 化 、 水 文 变化 等 方面 的 研究 提 
MERA (Promethei Lingula) 获取 的 探测 数据 ， 分 辨 4 


共 了 更 多 线索 。 文 [33] 利 ) 


2 埃 、 沙 占 比 导致 的 次 表层 介质 分 层 界面 


原 的 探测 揭示 了 层 状 沉 积 中 精细 的 分 层 结构 ， 为 火星 气 
jSHARAD 在 南极 高 原 普 罗 米 修 
8 南极 层 状 沉积 厚度 可 能 为 米 级 的 


围 为 1.7~$.3 us. 4.3-5.6 us, 5.6-7.1 us、11.0~12.9 


上 极 高 原 的 次 


几 个 次 表层 分 层 界面 (图 13)， 回 波 时 延 范 

us, 用 水 冰 的 介 电 常 数 实 部 3.5 将 时 延 转 换 为 深度 后 , 各 分 层 界面 深度 为 130~420 m、340~450 
m, 450~570 m. 890-1030 m。 文 [3 和 1 使 用 SHARAD 第 5192 轨 探测 数据 绘制 出 4 
表层 地 层 图 , 雷达 图 像 中 能 分 辨 出 不 同 成 分 的 冰 、4 

(图 14) ， 且 能 探测 到 层 状 沉积 下 方 的 基底 层 ， 此 外 ， 还 探测 到 一 些 深 达 数 百 


AR BOR E RS, 


这 些 螺旋 槽 将 浅 层 的 反射 回 波 分 为 数 段 , 同时 ， 对 北极 高 原 边缘 的 基底 层 露头 进行 探测 证 实 


了 格 明 娜 舌 状 高 原 
日 勒 里 


斯 带 塔 斯 


说 明 层 状 沉积 的 质量 负担} 


(Gemina Lingula) 下 不 存在 基底 层 ， 此 处 的 北极 


斯 平原 构造 CVastitas Borealis Formation) 上 。 


ARTA PAA fE PL 


BOA GS RAEI P RE 


在 极 厚 的 弹性 岩石 圈 这 一 假设 。 


探测 到 的 基底 较为 平坦 ， 


， 与 MARSIS 的 探测 结果 相符 ， 支 持 存 


图 13 SHARAD 第 2202 轨 探 测 数据 局 部 放大 图 。 图 中 


[33] 


显示 了 4 个 次 表层 反射 体 ， 这 4 个 反射 体 将 不 同 介 


Fig.13 Enlarged view of SHARAD data from orbit 2202. 4 subsurface reflectors separate different dielectric 


sequences [33] 
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图 14 SHARAD 第 5192 轨 探 测 数据 ， 已 经 过 时 延 -深度 转换 ， 使 用 的 次 表层 介质 介 电 常数 为 3 1。 


Fig.14 SHARAD data from orbit 5192. Range time delay has been converted into depth by assigning real 


permittivity of 3 below the detected ground surface Pl. 


2 ”国内 外 未 来 火星 探测 雷达 载荷 


2020 年 ,中 国美 国 、 阿 联 酋 以 及 一 些 民营 航天 企业 都 已 发 射 探测 器 开始 火星 探测 任务 ， 
而 欧洲 原 定 于 2020 年 发 射 的 火星 探测 器 由 于 疫情 等 原因 推迟 到 2022 年 的 火星 探测 窗口 。 我 国 、 
美国 的 探测 器 都 携带 了 由 近视 器 搭载 的 次 表层 探测 雷达 , 此 外 我 国 的 轨道 器 上 也 搭载 了 次 表 
层 探测 雷达 。NASA 的 MARS 2020 任 务 巡 视 器 “角力 号 ”搭载 的 雷达 为 火星 次 表层 实验 雷达 
成 像 仪 CThe Radar Imager for Mars subsurface eXperiment) , 简称 为 RIMFAX; 欧 空 局 ExoMars 
2022 计 划 的 “ 罗 莎 琳 德 .富兰克林 ”号 火星 车 搭载 的 雷达 实验 代号 “火星 水 冰 与 地 下 沉积 物探 
Mj” (Water Ice and Subsurface Deposit Observations on Mars) ， 简 称 为 WISDOM。 


2.1 EIA eg 

RIMFAX 主 要 目标 为 对 火星 次 表层 进行 成 像 ， 获取 火星 地 表 下 的 地 质 特征 , 并 提供 关于 
次 表层 成 分 的 信息 , 探 地 雷达 获取 的 剖面 图 像 能 够 揭示 层 积 岩 过 去 暴露 在 火星 表面 上 的 历史 ， 
比如 撞击 、 风 蚀 、 河 流 冲积 痕迹 等 等 5 。 

RIMEFAX 的 天 线 安装 在 巡视 器 的 尾部 偏 下 的 地 方 ， 其 天 线 为 超 宽 带 蝴 蝶 结 形 缝 除 天 线 ， 
使 用 时 将 在 接收 模式 与 发 射 模式 间 来 回 切换 。RIMEFAX 的 工作 频段 介 于 SHARAD、MARSIS 
与 欧 空 局 WISDOM 之 间 ，150-1200 MHz， 采 用 门限 调频 连续 波 〈Gated-FMCW ) ,进行 一 次 
全 频段 扫描 需要 的 时 间 为 1-20 ms, 具体 由 其 工作 模式 决定 ， 且 工作 频段 被 分 为 3 个 单独 的 扫 
描 频 段 ，150-300 MHz. 300-600 MHz 和 600-1200 MHz。RIMFAX 对 次 表层 地 貌 进行 成 像 的 
深度 超过 10 米 ， 垂 直 空 间 分 辨 率 〈 距 离 向 分 辨 率 ) 优 于 30 厘 米 ， 沿 火星 车 行进 路 线 每 10 厘 米 
进行 一 次 水 平 采样 〈 方 位 向 分 辨 率 ) P9, 

由 于 同时 接收 来 自 地 表 与 更 深 的 次 表层 回 波 是 不 可 能 的 , RIMFAX 划 分 了 3 种 工作 模式 : 
表面 模式 、 浅 层 模 式 和 滩 层 模式 。3 种 工作 模式 的 主要 区 别 在 于 接收 机 的 动态 范围 窗 ， 表 面 
模式 下 天 线 自身 反射 以 及 火星 表面 反射 回 波 在 窗 内 , 浅 层 模式 去 除了 天 线 自 身 反 射 , 动态 范 
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围 窗 将 表面 反射 以 及 浅 层 反射 包 括 在 内 ， 深 层 模式 将 天 线 自 身 反 射 、 表 面 反射 、 浅 层 反射 的 


回 波 都 置 于 动态 范围 窗 之 外 5 。 


RIMFAX 在 2015 年 4 月 、2016 年 4 月 分 别 进行 了 两 次 地 面 验证 实验 , 地 点 都 在 挪威 斯 瓦尔 
巴特 群岛 ， 该 区 域 为 冷 温 复合 冰川 ,表面 温 度 低 ， 部 分 冰川 底部 存在 较为 温暖 的 分 层 。2015 
年 数据 应 用 了 布莱克 曼 - 哈 里 斯 窗 函 数 , 进行 补 零 处 理 , 经 傅 里 时 反 变 换 后 得 到 成 像 ( 图 15 )， 


而 2016 年 在 此 基础 上 还 进行 了 去 背景 处 理 (图 16) 。 


图 15 (A) 2015 年 RIMFAX 外 场 试验 所 得 120 KT; B) A 图 中 双 


(A) 


(B) 


线 结构 放大 后 中 1。 


Fig.15 (A) The figure shows a 120 m long radar profile from 2015 RIMFAX field test. (B) The figure gives a 


zoomed in version of Fig. 14A P?l. 


图 16 2016 年 RIMFAX 外 场 试验 所 得 剖面 图 


meters and the depth scale is about 1.5 meters below the surface reflection ! 


moraine 


(表面 模式 ) ; 长 度 约 为 40m， 表 面 以 下 深度 约 为 1.5m™1。 
Fig.16 Radar profile using the Surface Model from 2016 RIMFAX field test. The length of the profile is around 40 
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欧 空 局 ExoMars 2022 计 划 “ 罗 莎 琳 德 . 富 兰 


克 林 ” 号 火星 车 搭载 的 WISDOM 雷 达 科 学 目标 


——— 
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AERE He Hb RES RD BTR SE, PRSE ACRI WR FEE SPI, FI Gk 
星 车 上 装备 的 钻头 进行 配合 , 雷达 可 以 沿 火 星 车 行驶 轨迹 在 钻探 前 获取 次 表层 信息 。 其 探测 
深度 约 3 米 ， 垂 直 向 的 空间 分 辨 率 可 达 厘 米 级 。WISDOM 可 以 进行 全 极 化 测量 ， 将 提供 关于 
三 维 地 质 结构 、 地 层 学 以 及 近 地 表 的 次 表层 演化 历史 的 关键 信息 , RB SL ERE DeC 
信息 ， 避 免 销 头 损坏 。 

WISDOM 是 一 种 特 高 频 步 进 频 率 雷 达 ， 频 段 500 MHz-3 GHz， 带 宽 2.5 GHz， 其 通过 天 
线 发 射 大 量 连续 波 信号 , 每 个 连续 波 对 应 一 个 频率 , 持续 时 间 为 一 个 步 长 At, 取 值 200 us", 
系统 的 带宽 由 发 射频 率 的 数量 Nyeq 以 及 发 射频 率 间 的 间隔 决定 ，WISDOM 的 Npregq 标 称 值 为 
1001，Nrreg 决 定 了 WISDOM 的 距离 分 辨 率 ， 而 频率 间隔 则 决定 了 雷达 的 最 大 模糊 距离 5” 。 
接收 到 数据 后 ， 将 对 数据 进行 希 尔 伯 特 变换 《使 信号 出 现 90" 相 移 ， 由 于 希 尔 伯 特 变换 只 适 
用 于 单 频率 信号 , 而 步 进 频率 信号 每 个 时 刻 只 发 射 一 种 频率 的 脉冲 , 因此 适用 希 尔 伯 特 变换 )， 
反 演 出 被 雷达 才 盖 环境 的 复 传 递 函 数 ， 再 进行 反 傅 里 叶 变 换 将 信号 转变 回 时 间 域 。 此 外 ,还 
将 对 回 波 信 号 施加 不 同 的 频 域 窗 函 数 以 减 小 距离 向 的 旁 渡 ， 进 行 滤波 以 提高 信号 的 信 吕 比 
(图 17) 。 
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图 17 从 左 到 右 顺序 .频率 域 测 得 信号 ， 加 窗 、 滤 波 减 小 旁 六 、 提 高 信 品 比 后 的 信号 ; (FFT 后 ) 时 间 域 合 
成 响应 ;时 间 域 合成 响应 (分 贝 表示 m1。 


Fig.17 From left to right: (a) Measured signal in the frequency domain. (b) Signal in the frequency domain after 


windowing and filtering to reduce side lobes and enhance the SNR. (c) Synthetic response in the time domain 


(after FFT) linear scale. (d) Synthetic response in the time domain in decibels "*l. 

随 着 巡视 器 行进 , 雷达 接收 到 的 数据 经 过 处 理 后 可 以 得 到 次 表层 地 质 结构 的 二 维 剖面 图 
此 时 的 数据 产品 只 提供 来 自 各 反射 体 回 波 的 时 延 。 经 过 对 巡视 器 行进 路 线 的 进一步 探测 ， 并 
使 用 一 个 金属 盘 〈 完 美 已 知 目标 》 进 行 定 标 ，WISDOM 可 以 获得 理想 情况 下 的 表面 回 波幅 
度 ， 并 估计 出 表层 的 介 电 常数 ， 确 定 表 层 的 厚度 。 

WISDOM 进 行 了 大 量 的 外 场地 面 验证 实验 ， 选 择 地 区 以 干燥 或 寒冷 环境 为 主 ， 包 括 智 


m 
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利 阿 塔 卡 马 沙漠 ， 意大利 埃 特 纳 火 山 ， 阿尔 捍 斯 山脉 冰川 以 及 澳大利亚 达 赫 施 泰 因 冰 穴 。 在 
多 种 环境 下 进行 实验 采集 到 的 数据 都 证 明 WISDOM 有 能力 分 辨 浅 层 精细 地 质 结 构 ， 为 研究 
地 理 、 地 质 演 化 历史 提供 线索 〈 图 18) 。 以 在 阿 塔 卡 马 沙漠 某 处 被 掩埋 的 河床 露头 进行 的 实 
验 数 据 为 例 ，WISDOM 成 功 探 测 到 河床 底部 以 下 约 50 cm 处 的 沉积 物 分 层 中 1。 


Distance (m) 
图 18 阿 塔 卡 马 沙漠 外 场 实验 WISDOM 雷达 图 像 ， 显 示 表 面 下 约 50 cm 深 处 为 河床 底部 沉积 物 的 界面 ， 被 


掩盖 的 河床 深 约 50 cm， 宽 约 7 mb 


Fig.18 WISDOM field test radargrams that reveal a small 50-cm-deep and 7-m-wide dried channel bed in the 


Atacama Desert (Chile) at about 30m from the left Pl. 


在 2010 年 10 月 意大利 埃 特 纳 火 山 进行 的 外 场 验 证 实验 开始 前 实验 场地 出 现 降水 , 造成 实 
验 场 中 火 成 碎 眉 沉 积 物 含水 量 较 高 ， 但 这 并 没有 严重 影响 到 WISDOM 发 挥 其 探测 能 力 。 
WISDOM 测 得 回流 数 据 绘 制 的 剖面 图 中 清晰 显示 出 次 表层 界面 的 形状 ， 齐 面 图 显示 的 最 大 
深度 为 3 m， 探 测 间隔 为 10 cmt 1。 由 于 WISDOM 天 线 的 方向 图 具有 方向 性 ， 可 进行 两 种 同 
极 化 探测 的 配置 有 助 于 WISDOM 绘 制 出 地 表 下 界面 的 三 维 形状 ， 确 定 其 所 在 位 置 范围 ， 埃 
特 纳 火山 以 及 达 赫 施 泰 因 冰 穴 等 地 的 实验 数据 处 理 后 的 产品 都 证 明了 这 一 点 (图 19) P". 


Depth (m) 


Distance (m) 


图 19 达 赫 施 泰 因 冰 穴 实验 数据 进行 三 维 重 建 得 到 的 基 岩 3D TR, Hea a TEL 1.5 m 厚 的 冰 层 下 5 。 


Fig.19 3D reconstruction of bedrock beneath a 1.5-m-thick ice layer (ice caves in Dachstein) P”. 
2.2 国内 雷达 载荷 
2020 年 7 月 23 日 发 射 的 我 国 火 星 探 测 计划 “天 问 一 号 ” 极 具 挑战 性 地 将 巡视 器 与 环绕 器 
一 起 送 往 火 星 , 巡视 器 上 搭载 了 中 国 科 学 院 电子 所 研制 的 火星 车 次 表层 探测 雷达 , 环绕 器 上 
搭载 了 中 国电 子 科技 集团 第 三 十 八 研究 所 研制 的 环绕 器 次 表层 探测 雷达 。 两 部 雷达 都 采用 了 
双 频 双 极 化 的 配置 , 同时 拥有 较 深 的 探测 深度 与 较 高 的 精度 , 双 极 化 探测 能 力 则 可 能 在 火星 
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次 表层 水 冰 探 测 方面 开拓 出 一 条 新 的 道路 。 
2.2.1 火星 车 次 表层 探测 雷达 


火星 车 次 表层 探测 雷达 (Rover Subsurface Penetrating Radar, RoSPR) 用 于 进行 火星 巡视 
区 地 表 、 次 表层 的 超 宽带 全 极 化 探测 ， 发 射 信号 为 线性 调频 信号 后， 高 频 模式 下 进行 全 极 
化 探测 。 低 频 通 道中 心 频率 55 MHz， 带 宽 80 MHz， 脉 宽 为 14$ us， 脉 冲 重 复 周 期 为 5 ms; 
高 频 通 道中 心 频率 为 1300 MHz， 带 宽 1700 MHz， 脉 宽 为 4 ms， 脉 冲 重 复 周期 为 5 ms。 以 雷 
达 工 作 参 数 推算 ， 在 低频 模式 下 ， 最 大 探测 深度 为 1125 m MERKE =4) ， 厚 度 分 辩 率 
为 0.94 m〔 介 电 常 数 : =4) ;在 高 频 模式 下 ， 最 大 探测 深度 为 30 m( 介 电 常 数 e =4)， 厚 度 分 


辩 率 为 4.4cm。 


低频 模式 使 用 的 天 线 为 两 组 单 极 子 天 线 , 安装 在 火星 车 顶板 下 方 。 高 频 模式 使 用 的 天 线 
为 Vivaldi 天 线 ， 安 装 在 火星 车 前 侧 板 上 ， 且 两 组 天 线 极 化 方式 为 相互 正 交 路。 高 频 、 低 频 
天 线 都 是 接收 天 线 、 发 射 天 线 分 开 ， 且 高 频 天 线 每 一 组 天 线 由 两 个 天 线 单元 组 成 ， 两 个 发 射 
天 线 单元 互相 正 交 ， 分 别 定 义 为 酉 极 化 发 射 天 线 与 V 极 化 发 射 天 线 ; 接收 天 线 同样 两 个 单元 
互相 正 交 安装 ， 定 义 为 H 极 化 接收 天 线 与 V 极 化 接收 天 线 。 因 此 ， 火 星 车 雷达 在 高 频次 表层 
探测 中 将 以 这 种 形式 进行 全 极 化 探测 ， 两 个 发 射 天 线 依 次 交替 发 射 正 交 极 化 的 电磁 波 ， 两 个 
接收 天 线 也 选 通 接 收 相 应 极 化 方式 的 回流， 可 构成 HH,HV,VH,VV 的 全 极 化 探测 (图 20〉。 


而 低频 通道 的 天 线 极 化 方式 为 HH， 无 法 进行 极 化 探测 。 


火星 车 次 表层 探测 雷达 进行 了 多 项 验证 试验 以 验证 其 工作 性 能 ， 并 研究 数据 处 理 方法 。 


实验 室内 场 试验 设置 了 一 个 长 7 m, 宽 3 m， 深 2.5 m 的 试验 池 ， 其 中 装 满 了 人 造 火 山 灰 〈 介 电 


常数 e = 2.94 土 0.56)， 并 在 深度 不 足 0.5 m 处 埋设 了 金属 管 、 塑 料 管 、 石 块 等 物体 (图 21 


) 。 


由 于 场地 限制 ， 该 试验 只 使 用 高 频 控 测 模式 ， 试 验 数据 去 除了 车 身 造 成 的 多 次 反射 ， 并 进行 
了 去 背景 操作 ， 数 据 处 理 结果 中 可 识别 出 金属 管 、 塑 料 管 与 池 底 的 回 波 《〈 图 22) UU. KE 


车 次 表层 探测 雷达 在 2018 年 、2019 年 前 往 陕 北 榆林 黄土 、 祁 连 山 老 虎 沟 冰 川 、 内 蒙古 乌 兰 哈 


达 火 山地 质 公园 等 地 进行 外 场 验证 试验 ,测试 其 高 频 、 低 频 模式 对 多 种 分 层 结构 的 探测 怕 


L Sk 
t HE o 


与 安装 在 “ 玉 免 ?系列 月 球 车 上 的 探 月 雷达 相 比 , 火星 车 次 表层 探测 雷达 不 再 使 用 脉冲 信和 号， 
而 改 为 线性 调频 信号 , 这 也 使 得 火星 车 次 表层 探测 雷达 的 数据 处 理 流程 与 探 月 雷达 有 所 不 同 。 


接收 天 线 (低频 ) 
Rx antenna(low frequency) 


发 射 天 线 (低频 ) 


Tx antenna(low frequency) 


接收 天 线 (高 频 ) 
< Rx antenna(high frequency) 
B RIKE QR) 


图 20 高 频 、 低 频 天 线 安装 示意 图 


Fig.20 The location of high frequency and low frequency antenna installation of RBMSPR. 
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图 21 填 满 人 造 火 山 灰 的 试验 池 中 1。 


Fig.21 the test pool filled with artificial volcanic ash", 


Time(ns) 


Metal pipe Plastic pipes 


Bottom of 
the pool 


Trace 


图 22 JU jd a A IE TRI PER AS D I T FL 


Fig.22 profile obtained when the rover on the pool top moving from the north to the south'!!. 

2.2.2 环绕 器 次 表层 探测 雷达 

我 国 搭载 在 “天 问 1 号 ”环绕 器 上 的 轨道 器 次 表层 探测 雷达 全 称 为 火星 环绕 器 次 表层 控 
测 雷 达 (Mars Orbiter Subsurface Investigation Radar, MOSIR ) . MOSIREMARSIS, SHARAD 
有 很 多 相似 的 地 方 ， 也 存在 很 大 差异 。MOSIR 采 用 了 双 极 化 探测 ， 这 种 探测 方式 可 以 通过 
计算 不 同 极 化 方式 回 波 的 斯 托 克 斯 参量 , 探测 水 冰 的 存在 , 还 可 能 提供 关于 火星 冰 盖 内 部 性 
质 等 的 线索 ,这 是 单 极 化 雷达 无 法 做 到 的 ,其 工作 模式 分 为 其 低频 探测 模式 、 低 频 探测 模式 、 
高 频 探测 模式 以 及 测 高 模式 ， 发 射 信号 都 为 线性 调频 信号 ， 各 工作 模式 中 心 频 率 、 带 宽 、 脉 
宽 、 接 收 时 窗 、 脉冲 重 复 频率 、 回 波 炙 加 数 以 及 采样 率 等 参数 不 同 , 见 表 2 所 示 。 与 MARSIS、 
SHARAD 相 似 的 是 ， 该 雷达 也 使 用 了 合成 孔径 技术 以 提高 方位 向 分 辨 率 。 

#2 三 种 模式 对 应 工作 参数 以 及 性 能 参数 


Table.2 Performance parameters and operating parameters of 3 modes of MOSIR. 


High Frequency Subsurface Sounding mode Low Frequency Subsurface Sounding mode Altitude Detection mode 


Pulse length(us) 60 70 (LF) /120 (HF) 120 


oT 
IV& ay ] 


ChinaX iv fERHT 
Data window duration(us) 95 120 (LF) /170 (HF) 170 
Center frequency (MHz) 40 12.5/17.5 40 
Frequency band (MHz) 20 5 20 
Sampling rate (MHz) 24 6 24 
Pulse repetition frequency 
810 741 1000 
(h=265km, Hz) 
Free-space range resolution 
12 48 
(m) 
Horizontal resolution in 
0.5-0.87 0.75-1.9 
along-track(km) 
Horizontal resolution in 
5-9 10-21 
cross-track(km) 
Penetration depth(km) 1.8 Ja 
Xi N ope EEL L— 4 Zu" NIE Ly Hn ez B i 
557 MOSIR 的 下 传 数据 给 数据 处 理 提 供 了 很 大 空间 ， 既 可 以 将 原始 数据 下 传 在 地 面 完成 肪 


558 ，” 压 、 聚 焦 等 处 理 ， 也 可 以 在 星 上 完成 成 像 后 下 传 。 以 原始 下 传 数据 为 例 ， 原 始 数据 下 传 且 完 
559 ”成 解 包 、 解 压缩 、 分 类 后 ,根据 雷达 地 面 大 回路 数据 以 及 地 面 定 标 数据 确定 收发 通道 、 天 线 
560 ”匹配 滤波 系数 ， 根 据 高 频 与 低频 模式 原始 回 波 数据 按照 文 [20] 给 出 的 方法 计算 出 低频 、 高 频 
561 ”相应 的 电离 层 补偿 系数 ， 点 乘 加 窗 后 的 标准 脉 压 系 数 ， 获 得 用 于 脉冲 压缩 的 参考 函数 ,一 次 
562 性 完成 脉 压 、 收 发 通道 补偿 、 天 线 响 应 补偿 、 电 离 层 补偿 。 完 成 脉 压 后 ， 按 照 目前 设计 的 数 
563 据 处 理 方法 ， 初 步 选 择 使 用 BP 算 法 完成 聚焦 成 像 处 理 。BP 算 法 实现 难度 较 低 ， 但 运算 量 较 
564 ”其 他 合成 孔径 处 理 算法 (如 CS 算法 等 ) 较 大 ， 其 原理 是 将 雷达 回 波 数据 投影 到 成 像 区 域 的 
565 ”每 个 像素 点 ， 每 个 像素 值 通 过 计算 相应 的 时 延 来 成 像 ， 具 体 实现 方法 为 在 每 个 FR 处 计算 环 
566 ” 绕 器 的 位 置 坐标 , 并 计算 21x300 个 聚焦 点 与 天 线 的 距离 , 以 距离 单元 (高 频 模 式 下 为 6.25m) 
567 ”为 单位 划分 每 个 聚焦 点 对 应 的 脉 压 后 数据 ， 以 波长 〈 高 频 为 7.3m， 即 光速 /中 心 频 率 ) 为 单 
568 ”位 划分 每 个 聚焦 点 对 应 的 相位 信息 ， 进 行 复 乘 、 多 普 勒 窗 函 数 处 理 后 进行 累加 处 理 ， 完 成 成 


569 fe 

570 MOSIR 除 了 利用 高 频 、 低 频 原始 数据 获得 对 应 的 电离 层 校正 系数 ， 还 可 以 利用 两 种 数 
571 ” 据 星 下 点 回 波 间 的 时 间 差 计算 出 对 应 的 电子 浓度 总 含量 TEC， 这 也 是 MOSIR 的 重要 数据 产 
572 品 之 一 。 

573 2019 年 ，MOSIR 在 内 蒙古 伴 口 县 进行 了 地 面 验证 试验 ， 该 实验 使 用 热气 球 巧 挂 环绕 器 


574 ”模型 以 及 安装 在 模型 上 的 次 表层 探测 雷达 进行 自由 球 飞 行 试验 以 及 系 留 飞行 试验 , 系 留 飞 行 
575 试验 主要 进行 天 线 隔 离 度 以 及 天 线 增 益 定 标 ， 自 由 球 飞 试验 使 用 测 高 模式 、 高 频 模式 、 低 频 
576 ”模式 进行 测量 ,验证 测 高 模式 测 高 精度 ， 经 地 形 矫正 、 时 间 校 正 后 测 高 精度 满足 10m 的 设计 
577 ”指标 要 求 。 高 频 模式 与 低频 模式 数据 用 于 数据 压缩 质量 验证 以 及 数据 处 理 流程 验证 , 结果 满 
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578 ”是 压缩 质量 要 求 ， 数 据 处 理 方法 与 流程 有 效 。 


579 
580 3 总 结 与 展望 

581 

582 由 于 雷达 在 较 低 的 频段 具有 穿 透 能 力 , 且 其 不 受 光 照 限制 , 可 以 对 天 体 不 受 光照 的 一 面 


583 ”或 永久 阴影 区 进行 探测 ， 在 对 月 球 、 行星 以 及 小 行星 的 探测 上 都 有 较 大 的 应 用 价值 ， 其 穿 透 
584 ”能 力 能 用 于 探测 星球 的 次 表层 结构 、 地 表 以 下 水 冰 的 存在 以 及 小 行星 的 内 部 结构 等 , 对 行星 
585 ”探测 研究 目标 天 体 的 地 质 、 气 修 、 水 文 演化 变迁 历史 具有 重大 意义 ， 且 雷达 可 具备 的 多 极 化 
586 ”探测 能 力 使 其 在 水 冰 以 及 干冰 探测 上 具有 独特 的 优势 。 但 由 于 雷达 通过 发 射电 磁 波 , 接收 物 
587 ” 体 反射 回来 的 电磁 波 以 进行 探测 , 本 身 的 性 质 决 定 其 需要 通过 一 系列 的 数据 处 理 才 能 反映 实 
588 ” 际 探 测 到 的 物理 量 , 日 进 行 轨 道 次 表层 雷达 探测 的 电磁 波 还 会 受到 行星 电离 层 的 影响 , 降低 
589 ”雷达 的 探测 精度 。 搭载 在 轨道 器 上 进行 次 表层 探测 的 雷达 一 般 采 用 合成 孔径 雷达 技术 , 需 
590 ”经 过 距离 向 和 方位 向 处 理 才 能 将 回 波 转化 为 可 读 的 图 像 ， 而 方位 向 处 理 方法 众多 ， 从 合成 孔 
591 ， 径 雷达 技术 起 步 至 今 , 已 出 现 chirp scaling 算 法 、 波 束 锐 化 算法 、ok 算 法 、RD (距离 -多 普 勒 ) 
592 ”算法 等 方法 , 目前 chirp scaling 算 法 与 波束 锐 化 算法 被 应 用 于 SHARAD 与 MARSIS 的 数据 处 理 
593 上， 而 相位 梯度 自 聚 焦 算 法 以 及 对 比 法 用 于 电离 层 校正 。 
594 MARSIS 与 SHARAD 是 目前 已 在 火星 轨道 上 工作 的 两 部 次 表层 探测 雷达 , 除了 通过 数据 
595 ”处 理 对 回 波 进行 成 像 , 在 后 期 的 研究 中 , 需要 使 用 MOLA 探 测 获得 的 火星 表面 数字 高 程 模型 
596 ”对 火星 表面 进行 建 模 ,用 于 火星 表面 杂 波 的 仿真 ; 使 用 多 种 方法 估计 火星 表面 以 及 次 表层 物 
597 ，” 质 的 介 电 和 常数、 损耗 角 正 切 等 数据 ， 以 此 为 依据 推测 地 表 之 下 的 物质 成 分 、 水 冰 的 含 侍 量 等 
598 ”科学 信息 。 
599 美国 于 2020 年 7 月 30 日 开始 火星 探测 任务 ， 其 巡视 器 携带 一 台 和 雷达。 我 国 首 次 火星 探测 
600 ”任务 在 2020 年 7 月 23 日 发 射 成 功 ， 包 括 环 绕 器 与 巡视 器 ， 各 搭载 一 部 雷达 。 我 国 的 环绕 器 次 
601 ”表层 探测 雷达 拥有 双 频 、 双 极 化 探测 的 能 力 ， 虽 然 是 我 国 首次 将 此 类 雷达 应 用 于 行星 探测 ， 
602 但 有 MARSIS 与 SHARAD 多 年 数据 处 理 的 技术 沉淀 以 及 其 对 火星 大 部 分 地 区 的 高 质量 探测 ， 
603 ”我 国 轨道 器 次 表层 探测 雷达 具有 在 MARSIS 与 SHARAD 的 基础 上 取得 更 多 成 就 的 潜力 。 
604 “MARSIS 与 SHARAD 对 火星 的 探测 覆盖 广 ,成 果 多 ， 如 何在 前 人 已 取得 如 此 成 就 的 基础 上 更 
605 ”进一步 是 我 国 雷达 研究 值得 思考 的 一 点 。 巡 视 器 雷达 虽然 更 改 了 体制 , 但 拥有 “玉兔 ”1 号 、 
606 ”2 号 的 技术 经 验 ， 以 及 月 球 车 雷达 的 优秀 表现 ， 我 们 应 该 对 火星 巡视 器 雷达 抱 有 信心 。 相 信 
607 ”我 国 首 次 火星 探测 任务 搭载 的 两 部 雷达 能 给 我 们 带 回 更 多 关于 火星 次 表层 结构 、 水 冰 等 的 科 
608 ”学 信息 ， 在 行星 科学 的 研究 上 发 挥 至 关 重 要 的 作用 。 
609 
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Overview of Mars Orbiter Subsurface Investigation Radar Data 


Processing Technology and Research Using Radar 
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Abstract: Mars is one of the most important targets of extraterrestrial exploration. 
Investigation and research on mars surface indicate that liquid water once existed on the mars 
surface, and water is basis and precondition of life. Therefore, one of today's mars exploration 


scientific objectives is searching for water both solid and liquid in mars subsurface. For the last 17 


years , the Mars Advanced Radar for Subsurface and Ionosphere Sounding (MARSIS) on ESA's 


Mars Express , and the Shallow Subsurface Radar(SHARAD) on NASA' s Mars Reconnaissance 


Orbiter has probed the mars subsurface at a global scale. Because of operating principle of radar, 
raw data of subsurface investigation radar should be processed in range and azimuth and correct 
ionospheric distortion in order to generate radargram. Many scientific achievements have been 
made in extracting scientific information from radargram. Data processing technology played a 
indispensable role in mars subsurface investigation, which shows great representativeness and 
reference value. We described the investigation and research about Martian subsurface with radar 
in this overview, introduced data processing technology using in Mars subsurface investigation, 
and presented part of findings of MARSIS and SHARAD. We also introduced several future 


Martian subsurface radars in this paper. 
Key words: Radar; Mars; Subsurface; range process; azimuth process; ionospheric 


distortion correction. 


